Exact diagonalization study of the tunable edge magnetism in graphene 
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The tunable magnetism at graphene edges with lengths of up to 48 unit cells is analyzed by an 
exact diagonalization technique. For this we use a generalized interacting one-dimensional model 
which can be tuned continuously from a limit describing graphene zigzag edge states with a ferro- 
magnetic phase, to a limit equivalent to a Hubbard chain, which does not allow ferromagnetism. 
This analysis sheds light onto the question why the edge states have a ferromagnetic ground state, 
while a usual one-dimensional metal does not. Essentially we find that there are two important 
features of edge states: (a) umklapp processes are completely forbidden for edge states; this allows 
a spin-polarized ground state, (b) the strong momentum dependence of the effective interaction 
vertex for edge states gives rise to a regime of partial spin-polarization and a second order phase 
transition between a standard paramagnetic Luttinger liquid and ferromagnetic Luttinger liquid. 

PACS numbers: 73.21.-b,75.70.Rf,81.05.ue,73.22.Pr 



I. INTRODUCTION 

Since it has first been isolated in the laboratory^ 
graphene, a two-dimensional honeycomb lattice of car- 
bon atomsj^ attracts much attention. In fact, graphene 
has multiple amazing properties. To name only a few 
of them, it ranges among the mechanically strongest 
materials,— it shows a quantum Hall effect at room 
temperature^ and, due to its unusual Dirac band struc- 
ture, it allows the study of relativistic quantum physics 
in a solid state environment^ Furthermore, its poten- 
tial application as the basis of the next generation of 
electronic devices stimulated great efforts to gain exper- 
imental control as well as theoretical understanding of 
this astonishing material. 

Usually a strong electron confinement increases the 
strength of electron-electron interactions by pushing the 
electrons close together. However, in spite of the ex- 
treme electron confinement to only one single layer of 
atoms, many experiments in graphene may be explained 
by assuming the electrons to be non-interacting. This is 
especially true for experiments probing the bulk proper- 
ties of graphene, as the bulk density of states vanishes at 
the Fermi level, suppressing the manifestation of interac- 
tion effects. On the other hand, the properties of zigzag 
edges differ greatly from the bulk properties. So-called 
edge states, i.e. one-dimensional states with very small 
bandwidth, localized at these edges, give rise to a peak 
in the local density of states at the Fermi energy^ The 
enhanced density of states allows the electron-electron 
interaction to drive the zigzag edges to a ferromagnetic 
state with a magnetic moment localized at the edge^ - — 
This phenomenon is known as edge magnetism. 

At normal graphene edges the electron-electron inter- 
action is so strong and the bandwidth of the edge states 
is so small that the spins of all electrons in the edge states 



are completely aligned. However, as has been proposed 
recently, graphene / graphane interfaces provide means to 
tune the bandwidth of the edge states to regimes in which 
the edge starts to depolarize and the edge magnetism is 
gradually suppressed until, for a critical edge state band- 
width, the magnetism disappears^ 

In Ref. [TH it was argued that this interaction-induced 
magnetism can be understood on the basis of an effective 
model, describing the interacting one-dimensional edge 
states only, while the bulk states are neglected. What 
at first glance appears to be a contradiction to the Lieb- 
Mattis theorem^ stating that the ground state of inter- 
acting electrons in one dimension cannot be spin polar- 
ized, can be resolved by noting that the effective edge 
state model does not fulfill the prerequisites of the Lieb- 
Mattis theorem The deeper reasons for the existence of 
a ferromagnetic ground state in a one-dimensional inter- 
acting electron system, however, remained elusive. The 
present work is devoted to this issue. 

In this paper we present a systematic exact diagonal- 
ization analysis of interacting edge states. Two striking 
features of edge states turn out to be most important 
for their magnetic properties: (a) the edge states exist 
only in a restricted part of the Brillouin zone and (b) the 
transverse edge state wave function has a strong char- 
acteristic momentum dependence. These features have 
consequences for the effective low-energy theory, namely 
(a) no umklapp processes are allowed in the interaction 
Hamiltonian and (b) the interaction vertex acquires an 
unusually strong momentum dependence. In order to be 
able to study the consequences of these two features, we 
introduce a generalized model in which we add an arti- 
ficial interaction term describing umklapp processes and 
allow the momentum-dependence of the interaction ver- 
tex to be tuned from a momentum-independent vertex, 
as in usual metals, to the full momentum-dependence, 
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as it is found in edge states. Therefore, the generalized 
model can be tuned continuously from the limit in which 
it describes edge states to a limit which corresponds to 
usual one-dimensional metals such as the Hubbard chain. 
We solve this generalized model for graphene zigzag edges 
of finite length L = 48 unit cells (i.e. ~12 nm) by exact 
diagonalization using the Lanczos method for the deter- 
mination of the ground state of the effective modelJ^— 
The paper is organized as follows. In Sec. [IT] we review 
the direct model as it has been derived in Ref. UH and 
introduce the more versatile generalized model with ad- 
ditional tunability. In Sec. IIIII the exact diagonalization 
analysis of the generalized model is presented. Finally, 
the results are discussed in Sec. Uvl 



II. EDGE STATE MODELS 

In this section we introduce the models on which our 
analysis is based. The edge state model obtained from 
the direct projection of the honeycomb lattice Hubbard 
model to the Fock space spanned only by the edge states 
has been discussed in Ref. This model will be called 
the direct model in the following. We identify two impor- 
tant features of the direct model: (a) the restriction of 
the Brillouin zone for the edge states and (b) the strong 
dependence of the transverse localization length on the 
momentum along the edge. After having analyzed the 
consequences of these features for the effective interac- 
tion vertex, we propose a generalized edge state model in 
which these features can be tuned. This allows us to in- 
vestigate the impact of each of these edge state features 
on the magnetic properties. In particular, the general- 
ized model can be tuned continuously from a Hubbard 
chain limit, i.e., a usual one-dimensional metal without 
any ferromagnetic ground state, to the edge state limit 
with its ferromagnetic ground state. 



A. Direct derivation from the honeycomb model 

We start from the simplest possible non-interacting 
tight-binding model of electrons in graphene zigzag rib- 
bons, taking into account only nearest neighbor hoppings 
of 7r electrons % = Y^Uj) a c la c j<^^ wnere (hj) runs over 
nearest neighbor sites of a half-infinite honeycomb lattice, 
i = (m, n, s) is a collective site index for the (m, n)th unit 
cell and the s = A,B sublattice (see Fig. [T]), and Ci a an- 
nihilates an electron at site i with spin a. Since we are 
exclusively interested in the zero energy eigenstates of H, 
the actual energy scale of % is unimportant so that we 
may drop itJ^ The zero energy states are created by the 
fcrmionic edge state operator 
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where u p = — 1 — e lp , p is the momentum in m direction 



with the number of unit cells L mm direction, i.e., along 
the edge. The p-dependent normalization constant M p = 
\J 2 cos(p — 7r) — 1 can be interpreted as the weight of the 
edge state wave function right at the edge atoms where 
n = 0. It is easily seen that He pa = 0. As the edge state 
wave function is only non-zero on the B sublattice we 
omit the sublattice index, setting it to s = B. 




FIG. 1. (Color online) Part (a) shows the half-infinite hon- 
eycomb lattice. The solid elipses (green) are the complete 
unit cells in the bulk region and the dashed elipses indicate 
the cut unit cells at the a edge. The n and m directions 
are indicated as well as the sublattice indices A, B. Part (b) 
shows the modulus square of the transverse (in n direction) 
edge state wave function \ip p (n)\ 2 on the B sublattice sites in 
the reduced Brillouin zone ^ < p < 4^. The extreme mo- 
mentum dependence of the localization of %p p in n direction is 
crucial for the magnetic properties of edge states. 

The two most important features of the edge state wave 
function [Eq. (UJ)] are : (a) The edge state only exists for 



In the rest of the Brillouin 



(along the edge), and 



-*pn<7 



-ipm 



momenta ^ < p < 
zone the edge state wave function is not normalizablc, 
as \u p \ > 1 for these momenta, (b) In n direction the 
edge state is sharply localized at the edge for p = w, 
wereas for p close to one of the Dirac points K = 
and K' = ^, the wave function delocalizes into the bulk 
[see also Fig. H{b)]. These two edge state properties 
are stable against adding more details, such as second- 
nearest neighbor hopping or various edge passivations, to 
the honeycomb Hamiltonian H i 13 ' 20 The detailed analy- 
sis presented in this paper will clarify that the existence of 
edge magnetism and in particular its tunability are con- 
sequences of these two edge state properties. 

The p-dependence of the localization length of if> p has 
consequences for the edge states' self-energy e as well 
as for their interaction vertex function T. Neglecting 
the bulk state contributions^ the self-energy correction 
due to a perturbation V , which is invariant along the 
edge, is given by eo{p) = (ip p \V\tp p ). Due to the delo- 
calization of ippin) for p near K,K', edge- localized per- 
turbations V lead to self-energy corrections for which 
e (K) = eo(K') = while e (n) - ||V||. For suffi- 
ciently well behaved perturbations the self-energy correc- 
tion gives rise to a smooth edge state energy dispersion 
with a bandwidth A ~ ||V"||. For a large class of these 
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edge-localized perturbations, the self-energy correction 
approximately has the form 

e (p)^^A, (2) 

with Afp the p-dependent weight of the edge state wave 
function right at the edge. Eq. ^ expresses that an 
edge state which is more localized at the edge experiences 
a stronger self-energy correction from an edge-localized 
perturbation than an edge state which is dclocalized into 
the bulk region. Examples of such perturbations are 
edge passivations, graphane termination^ or local in- 
teractions with the substrate. Note that the edge state 
bandwidth A is experimentally tunable in various ways 
so that we consider A as a free parameter. Therefore, the 
nonintcracting part of the direct model (dm) edge state 
Hamiltonian is given by: 

tf dm = -A]T£X 2 4,^- ( 3 ) 

a p 

where the sum is restricted such that only edge state 
operators e pa with ^ < p < ^- appear. 

The effective interaction of the edge states, derived 
by projecting the Hubbard Hamiltonian on the two- 
dimensional honeycomb lattice Hu = Uj^i c !t c it c lj. c j| 
to the Fock space spanned by the edge states, readsi^ 

Hfn = |^'r(p ) p', g )et +3t e pt e+,_ g ^ i . (4) 

Together, we have the effective Hamiltonian H = 
Hq m + Hf m . The interaction vertex is given by the over- 
lap of the wave functions of all four fermions (with mo- 
menta p + q,p-,p' — q,p') participating in the interaction 

oo 

F(p,p',q) = ^%+ q {n)i>p{n)^l<- q {n)'<l ! p'{n) 

= N p+q N p N pl - q N p > . 

1 — U p+q U P U p' -q U P' 

While the denominator, resulting from the geometric se- 
ries over n, turns out to lead only to unimportant quanti- 
tative corrections, the numerator of T, which is the prod- 
uct of the wave function weights at the edge N p for each 
of the four fermion operators, leads to the momentum- 
dependence of the interaction strength which is impor- 
tant for the stability of the weak edge magnetism. Es- 
sentially, the effective interaction becomes stronger the 
more localized the participating fermions are, i.e., the 
closer their momenta are to p = it. If one or more of the 
momenta are close to the Dirac points where the 

edge state wave functions delocalize into the bulk, the 
effective interaction is suppressed (see Fig. [J). Note that 
setting the denominator in Eq. ([3]) to unity corresponds 
to assuming that the Hubbard interaction is only present 
at the outermost line of carbon atoms right at the edge. 



Such an approximation has been used in Ref. [U We 
find that this approximation is inessential for the edge 
magnetism, leading only to quantitative corrections. 

An important consequence of the restriction of the p 
summation in Eq. (j4]) is the absence of umklapp pro- 
cesses. As explained above, edge states only exist in one 
third of the Brillouin zone, i.e. for ^ < p < so 
that no four fermion process with momentum ±2tt ex- 
ists. Indeed, within the restricted Brillouin zone, the 
process with the largest possible total momentum p iot 

is e L/3T e27r /3T e L/3i e2 V3i> i- c - Ptot = f < 2tt. Pro- 
cesses with larger total momentum leave the restricted 
Brillouin zone and are therefore suppressed, as they in- 
volve the overlap of edge states and bulk states, which 
is small. Also, most of the bulk states live in a different 
energy regime than the edge states^ 

Thus, we have identified two properties of edge states 
which make them fundamentally different from usual one- 
dimensional conductors: 

(a) Due to the restricted Brillouin zone, umklapp pro- 
cesses are forbidden. 

(b) The transverse localization Af p of the edge state 
wave function gives rise to a tunable band width 
eo{p) — A/'pA. Furthermore, the interaction vertex 
becomes weaker if the momenta of the participating 
fermions approach a Dirac point, i.e. T(p,p',q) cx 

Afp+qAfpAfp'-qAfp'. 

We will show that these properties are the basis for the 
magnetism at graphene edges. 

B. Generalized model 

We now introduce a generalized model in which the 
different aspects of the effective electron-electron inter- 
action, found in the previous subsection, may be tuned 
independently. For this, we map the edge state opera- 
tors e pM which correspond to right (left) moving modes 
for 7r < p < 4^ (^ < p < 7r), to fcrmionic operators 
c-kro in which r = R, L specifies the direction of mo- 
tion and -f < & < f , i.e. c kra = e fc+7r+r7r/6iCT and 
p = k + 7r + m/6. The direction of motion r = R, L cor- 
responds to r = ±1 when used in formulas. Note that the 
zero point of k has been shifted so that k = corresponds 
to p = 7r±7r/6 for right and left movers, respectively (sec 

Fig. HD. 

For the non-interacting part of the generalized edge 
state model we assume a linear spectrum with slope ±vf 

tt/6 

#0 = V F ]T ^ ( rk ) C kra C kra- (6) 
r=R,L fc=-7r/6 

<T=t,4- 

This linearization of the self-energy [Eq. @] only leads 
to inessential quantitative corrections (see also Appendix 
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FIG. 2. (Color online) Left: The self-energy ~ A/" p 2 of the 
direct model (dashed line) and the linearized self-energy (solid 
lines). Right: The four possible interaction processes. The g$ 
process is not allowed for graphene edge states. 



The partitioning into left- and right-movers always 
gives rise to four terms in the interaction part H\ of 
the Hamiltonian involving different combinations of left- 
and right-moving modes. Conventionally, these terms 
are called <7i, <?2 ; 9z, 9 a processes (see Fig. [2] and Ref. |22| ). 
92,9a correspond to forward scattering, involving pro- 
cesses that scatter only between modes with the same di- 
rection of motion, g\ refers to backward scattering, and 
<?3 are the umklapp terms which are forbidden in edge 
states. Note that, unlike in usual g-ology,— we may not 
assume that the coupling constants for the individual gi 
processes are constant. The momentum-dependence of 
the gi must be taken into account. 

The two forward scattering processes <?2 and 9a may 
be merged together into one Hamiltonian^ 



rrfs \ \ ^' nr or nr or' .> 

1 C fc+grt Cfer t C !'- (? r'; C fc'>-'l : > CO 



r,r' k,k' ,q 



where : A : enforces the normal order— of the operator A. 
The primed sum is restricted such that |fc| < 7r/6 for all 
momentum arguments k in the electron operators. In or- 
der to be able to change the amplitude of the momentum- 
dependence of the interaction vertex T(p,p',q) [see Eq. 
(O], we introduce the factors 



(8) 



from which we build the interaction vertex for the gen- 
eralized model. The factor Ti € [0,6/7r] quantifies the 
momentum-dependence. For T\ = the interaction is 
momentum-independent. This limit corresponds to usual 
one-dimensional Hubbard chains. For T± = 6/n the in- 
teraction goes to zero if at least one of the fcrmions is 
close to the upper band edge (k = m/6). This cor- 
responds to the direct model, where the trigonometric 
term under the square root in Af p has been replaced by a 
linear approximation. The differences between the gener- 
alized model in the edge state limit and the direct model 
only lead to quantitative rcnormalizations of the critical 
point, as shown in Appendix IBl The essential property 
of the interaction vertex is that it approaches zero if one 



of the fermion momenta gets close to the Dirac points. 
This feature is present in the direct and in the generalized 
model with T± = 6/n. 

The form of the backscattering (gi) Hamiltonian H^ s 
is similar to H[ s . However the scattering takes place be- 
tween left- and right-movers 



H ] > s = A 



U 



bs" 



EE 

r k,k' ,q 



or q—- 
D k+q D k 



-qSl'X 



C \+q,r^ C k,-rA C \' -q-r.\. C k' ,r,±- (9) 

We have introduced the additional parameter Ab s which 
allows us to tune the overall strength of the gi processes 
relative to the 92,9a processes. Abs = 1 corresponds to 
the physical backscattering strength which is required by 
SU(2) invariancei22 Nevertheless, we will investigate the 
consequences of a suppression of backscattering since this 
will be important for a bosonization analysis of the gen- 
eralized model which will be presented in an upcoming 
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papers 

As already pointed out, an important feature of edge 
states is the absence of umklapp processes in the effec- 
tive electron-electron interaction. However, in order to 
be able to compare the edge state model to a Hubbard 
chain, we add an artificial umklapp process with relative 
strength A us to the Hamiltonian of the generalized model 



A, 



\ \ qt Q—rQr q — r 



k,k' 



-k+q,r^ c k-r,t c k'-q,r,i, c k'-r,i- (10) 

Varying A us between 1 (Hubbard chain limit) and (edge 
state limit) allows us to investigate the consequences of 
the presence of umklapp processes for one-dimensional 
fcrromagnctism. 

Alltogethcr, the four parameters vp/U, Abs, A us , and 
Ti define the phase space of the generalized model 

H = H + H{ S + H^ s + H? s . (11) 

The following limits of this model may be identified: 

1. Edge state limit: the generalized model with the 
parameters Abs = 1, A us = 0, and Ti = 6/71", is a 
good approximation of the direct model. 

2. Hubbard chain limit: for Ab s = 1, A us = 1, and 
Ti = 0, the generalized model essentially describes 
a one-dimensional Hubbard chain. The only differ- 
ence is the assumption of a linearized single-particle 
spectrum instead of the 2 cos(fc) dispersion. 

Note that it is important to work in the k-space for- 
mulation because it is difficult to control the umklapp 
scattering or the momentum dependence of the interac- 
tion vertex in a real space formulation. One reason for 
this is that an interaction vertex T(p,p',q) with a non- 
trivial p, p' dependence does not transform to a real space 
interaction of the form V(x — x') but to a complicated 
non-local interaction. This also hampers the application 
of DMRG methods to this problem. 
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III. EXACT DIAGONALIZATION 

The ground state of the generalized model is calcu- 
lated for finite sized zigzag edges up to L = 48 by 
the Lanczos exact diagonalization metho d 2515 " —. The 
magnetic properties of the ground state depend on 
the ratio between the kinetic energy and the poten- 
tial energy vf^/U, which is experimentally tunable at 
graphene/graphane interfaces^ Three additional tuning 
parameters ri,A us ,Ab s , which are not accessible experi- 
mentally, have been added in order to be able to study 
the significance of the momentum dependence of the in- 
teraction vertex (Ti), the influence of the absence of umk- 
lapp scattering (A us ) in edge states, and also the impor- 
tance of backscattering (Abs)- With those artificial pa- 
rameters, the generalized model may be tuned continu- 
ously from a Hubbard chain limit to the edge state limit. 
In both limits the model describes an interacting one- 
dimensional metal. The magnetic properties in these two 
limits, however, differ strongly: while the usual Hubbard 
chain (with umklapp scattering and without momentum 
dependence) does not give rise to a ferromagnetic ground 
state, the edge states (without umklapp scattering and 
with momentum-dependent interactions) show two mag- 
netic phases in addition to the non-magnetic Luttinger 
liquid phase: for strong interactions the saturated edge 
magnetism2r— is recovered, while for intermediate inter- 
action strengths, a ferromagnetic Luttinger liquid ap- 
pears. 

The Hamiltonian H [Eq. (|Ilj) ] conserves the numbers 
iV-fjiVj. of up-spin and down-spin electrons, so that H is 
block diagonal in the S z subspaces, which we define by 
the total spin-polarization in z direction 

S 2 = i(iV t -A^) = 0,I,2,...,iV/2. (12) 

The total number of electrons N = Nf+N± = L/3 is kept 
constant. This corresponds to half-filling. Note, however, 
that the filling is physically relevant only if umklapp scat- 
tering is present (i.e. A us > 0). For the edge states in 
which we are finally interested, umklapp scattering is for- 
bidden so that the filling is irrelevant as it only leads to 
quantitative renormalizations of the interaction strength 
and the Fermi velocity. In the following, we determine 
the ground state of H in each S z subspace separately. 

Note that by the definition of the S z subspaces we have 
chosen a spin quantization axis. The Hamiltonian H, 
however, is SU(2) symmetric if the backscattering is at 
its physical strength Abs = 1- Furthermore, since we are 
dealing with finite systems, there will be no spontaneous 
rotational symmetry breaking. Thus, edge magnetism 
will become manifest in a (2S + l)-fold ground state de- 
generacy, corresponding to a high spin (S) state. For 
instance, if in a system with N = 2 electrons the low- 
est energy states in the subspaces S z = —1,0,1 are the 
degenerate ground states, then the h spins of two elec- 
trons point into the same direction, building an S = 1 
super spin. Because the SU(2) symmetry of the individ- 



ual electron spins is not broken, also this composite super 
spin has full rotational symmetry. The S z quantum num- 
bers of the degenerate S z subspaces then correspond to 
the magnetization of this composite spin system. Note 
that the spin-orbit interaction lowers the symmetry of 
the super spin, as it breaks the SU(2) invariance of the 
individual electron spins which form the super spin. 

For practical reasons, we extract the total spin quan- 
tum number S of the ground state from its ground state 
degeneracy (25 + 1), which is obtained from the S z sub- 
space ground state energies. We have checked that this 
is equivalent to calculating the total spin S of the ground 
state directly. 



A. Hubbard chain vs. edge states 

First we study the crossover from a usual Hubbard 
chain to interacting edge states. As explained above, 
the generalized model can be tuned continuously between 
these two limiting cases by means of the parameters A us 
and Ti. We postpone the analysis of backscattering to 
the next subsection and set Abs = 1 here. 

It is most instructive to begin with the Hubbard chain 
limit of the generalized model, which is characterized by 
the full umklapp process strength A us = 1 and a sup- 
pressed momentum dependence of the interaction Ti = 0. 
With this parameter set the direct model resembles a 
one-dimensional metal with a linear single particle dis- 
persion instead of a cos-dispersion^ The lowest eigen 
energies in the different S z subspaces for the parameter 
set described above are shown in Fig. [3j Obviously, the 
ground state is non-degenerate and resides in the S z = 
subspace for arbitrary vftt/U , except for the limit of in- 
finitely large U. Thus, as expected, no ferromagnetic 
phase transition exists for the Hubbard chain limit of the 
generalized model at finite vftt/U , in consistence with 
the Lieb-Mattis theorem^ which states that the ground 
state of a system of one-dimensional interacting electrons 
has zero total spin and is non-degenerate with higher spin 
subspaces ^2 

Next, the generalized model is tuned away from the 
Hubbard chain limit by suppressing the umklapp scat- 
tering A us < 1- Suppressed umklapp scattering is one 
of the properties of edge states which makes them fun- 
damentally different from usual one-dimensional metals. 
In the inset of Fig. [3J the lowest eigen energies of the 
S z = and the S z = ±7V/2 (full spin-polarization) sub- 
spaces are shown as A us is reduced from 1 to in steps of 
0.2. For any A us < 1 there is a nonzero critical value for 
vftt/U below which the lowest energy states of these two 
subspaces and also for all S z in between (not shown in 
the inset of Fig. [3j are equal. This corresponds to a high 
spin state of size S — 6 in which the spins of all electrons 
point into the same direction. The critical point at which 
the transition between 5 = and S = N/2 takes place 



6 



Lowest eigen energies in different S z subspaces 
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FIG. 3. (Color online) Lowest eigen energies in different S z 
subspaces for N = 12 in the Hubbard limit with Abs = 1, 
A us = 1, Ti = 0, and (7 = 1. The inset shows the lowest 
eigen energies of the S z — 0, 6 subspaces (zero and full spin- 
polarization) as the umklapp scattering is suppressed. The 
lower lines correspond to A us = 1 and the higher lines to 
Aus = 0. For the lines in between, A us decreases in steps of 
0.2. 



depends on the umklapp scattering strength 



VpTT 



u 



(r 1= o) 



oc(l-A us f. 



(13) 



For N = 12 we find for the exponent y ~ 0.5 ± 0.02. Ob- 
viously, the absence of umklapp scattering allows a high 
spin ground state. However, for A us < 1 and Ti =0, 
the system instantly jumps from zero polarization 5 = 
to the maximal possible polarization S = N/2 at the 

critical point [vftv/U]^ . This is a first order phase 
transition. For the case of completely suppressed umk- 
lapp scattering A us = this is shown in Fig. H](a), where 
the lowest eigen energies of all subspaces are plotted as 
a function of Vftt/U: the S z = subspace contains the 
non-degenerate ground state until at the critical point the 
lowest energy eigenstates of all subspaces form the de- 
generate ground state; no intermediate regime of vpir/U 
exists in which there is only a degeneracy between some 
of the S z subspaces. 

The reason for the instant jump in the total spin 
is as follows: once the Stoner criterion vpir/U > 

[vpiT /t/] c^it _ °' > ^ S me t, the interaction energy gain 5Eu(S) 
associated with developing a certain spin polarization S 
is larger than the corresponding kinetic energy penalty 
5Eki n (S). Unlike in two or three dimensions, however, for 
one-dimensional systems with momentum-independent 
interactions, SE(S) = SEu(S) + SE^^S) has no min- 
imum, i.e. SE(S + 1) < SE(S), for all S < S max . Thus, 
the system instantly 'flows' to the highest possible po- 
larization 5 max , once the Stoner criterion is met. This 
is a rather common feature of one-dimensional systems 
with a constant interaction vertex (such as the Hubbard 



FIG. 4. (Color online) Lowest eigen energies with com- 
pletely suppressed umklapp scattering A us = and differ- 
ent momentum-dependencies Ti. Furthermore, N = 16, 
Abs = 1, and U = 1. Part (a) shows the case of a momentum- 
independent interaction vertex (ri = 0), where the ground 
state degeneracy jumps from 17 (= 25 max + 1) directly to 
1 at the position indicated by the arrow. Part (b) shows 
the case of an interaction vertex with maximal momentum- 
dependence Ti = 6/7T. The arrows indicate a change in the 
ground state degeneracy. 



interaction) and can easily be observed in a variational 
calculation of the ground state properties (see Appendix 

The momentum-dependence of the interaction vertex 
(Ti > 0) reduces the interaction energy gain as the spin- 
polarization S becomes larger. This is because for larger 
S, the Fermi level of the spin-up right-movers is shifted to 
higher momenta where the interaction is suppressed by 
the Sl =R factors [see Eq. (J7J]. Similarly, for the spin- up 
left-movers, the Fermi level is then shifted to smaller mo- 
menta, where the Sl =L suppress the interaction.— As a 
result, SE(S) develops ci minimum at S — 'S'min < s 

max ' 

and the system is stable there. Intuitively, this may be 
understood on the basis of a variational calculation (see 
Appendix [A"|l . Within exact diagonalization one finds 
that with Ti = 6 /it, the ground state degeneracy in- 
creases successively from to 2S max + 1 by first adding 
the lowest energy eigenstates of the S z = ±1 subspaces 
to the ground space, and then the S z = ±2 subspaces 
and so forth. This is shown in Fig. |4] (b) . 

If the total spin S is plotted as a function of vpw/U, 
S decreases from to in steps. These steps corre- 

spond to the positions vftt/U, where the degree of the 
ground state degeneracy changes, indicated by arrows in 
Fig. |H For Ti = 0, there is one step where the spin- 
polarization jumps from S = SVax = N/2 to S = 0, 
while for the maximal Ti = 6/V, there are N/4 steps 
at each of which the spin-polarization is decreased by 
AS = 2.— Fig. [5] shows these two limiting cases and 
how the steps evolve as Ti is varied from to 6/V. The 
momentum-dependence must have a minimum strength 
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FIG. 5. (Color online) Dependence of the spin-polarization 
S on vftv/U for different strengths of the momentum- 
dependence Ti from Ti = (rightmost curve) to Ti = — 
(leftmost curve) in steps of AFi = y|^. The length of the 
edge is L — 48, Ab s = 1 and A us = 0. The smooth curve is a 
power law fit to the plateau centers of the Ti = — steps with 
exponent j3 = 0.5 (see text). 



Ti > T1 ~ 1, in order to break the one big spin- polariza- 
tion step of height N/2 into many small steps of height 
2. Thus, for T± > there is a regime of weak edge 
magnetism, meaning that the total spin 5* of the ground 
state is smaller than the maximal spin SVax, in addition 
to the usual saturated edge magnetism for small vpir/U 
(i.e. S = Smax) and the Luttinger liquid regime for large 
VFit/U with 5* = 0. Figure [6] shows a diagram in which 
the phase boundaries between the Luttinger liquid (LL), 
the saturated edge magnetism (SEM) and the novel weak 
edge magnetism (WEM) are shown for different system 
sizes N = 8, 12, 16. 

Note that the non-zero Tf found in the exact diag- 
onalization reveals a weakness of the fermionic mean- 
field theory in which this minimum momentum depen- 
dence, above which a WEM regime appears, is zero (see 
Appendix |A|. A non-zero rf means that the small 
momentum-dependencies which always follow from a de- 
pendence of the Bloch wave functions in usual one- 
dimensional conductors on the momentum are not neces- 
sarily sufficient to stabilize the weak edge magnetism; the 
momentum-dependence of the interaction vertex must be 
sufficiently strong for this. 

In the limit L — > oo, which cannot be accessed 
within exact diagonalization, of course, S/S max becomes 
a smooth function of Vpir/U. We approximate this 
smooth function by a power law 



S / S n 



Vftt' 

U / crit. 



IT 



(14) 



Fermionic mean-field theory (see Appendix [3} predicts 
f3 = 0.5. Because the exact diagonalization study is lim- 
ited to small systems N < 16, it is difficult to obtain a 



decent estimate of the exact exponent /3 within this work. 
Fitting the edge state limit of the generalized model to 
the center of the plateaus of the N = 16 results of the 
exact diagonalization, we obtain 0.44 < j3 < 0.61, depen- 
dent on how many plateaus are included in the fit. The 
critical point {^ff-) crit — 0.760 for this fit is obtained 
by extrapolating the rightmost step from the data sets 
TV = 4, 8, 12, 16 to7V = o. 



Phase diagram for A us = 
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FIG. 6. (Color online) Phase diagram for lengths L = 24, 
L = 36 and L = 48 . For small velocity dependence Ti of 
the interaction, only the Luttinger liquid (LL) phase and the 
saturated edge magnetism (SEM) phase exist, whereas above 
the critical value of Fi = Ff the weak edge magnetism (WEM) 
phase appears. 

Interestingly, Ti not only affects the order of the transi- 
tion but also the critical vf^/U. This is also not correctly 
predicted by the mean-field approach (see Appendix [A"| , 
which, independently of F\, finds vpTt/U — 1 to be the 
critical point. For small Y\ and N = 16 the exact diago- 
nalization gives 



u 



0.84- 0.17T?. 



(15) 



For the maximal T\ = 6/n, the position of the left- 
most step can be calculated by exact diagonalization for 
very large systems^ We performed calculations for sys- 
tem sizes up to I = 180 in order to extrapolate this 
step position. Within the limits of the accuracy of this 
extrapolation, the critical point vf^/U = 0.5 ±0.001 be- 
tween the SEM and the WEM regime coincides with the 
mean- field prediction (see Appendix [X]) . This extrapo- 
lation to the thermodynamic limit, in combination with 
the extrapolation of the critical point between the WEM 
and the LL regime, is a strong evidence for the existence 
of the WEM phase for 0.5 < v F n/U < 0.760 in the ther- 
modynamic limit. 

For completeness we note that our exact diagonaliza- 
tion analysis shows that a SEM phase also exists in the 
general model with umklapp scattering A us = 1 if Ti > 0. 
However even for Ti — 6/ir there is no weak edge mag- 
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nctism phase between the Luttinger liquid and the satu- 
rated edge magnetism as long as A us = 1. 



B. The relevance of backscattering 

The backscattering Hamiltonian H^ s is important for 
the SU(2) invariance of the Hamiltonian. It is easily 
seen that only for Ab s = 1 the SU(2) symmetry is pre- 
served. At real graphene edges, of course, the backscat- 
tering cannot be tuned experimentally. Nevertheless it 
is interesting to study the consequences of a suppression 
of H^ s since in a bosonization treatment of the general- 
ized model translates to a sine-Gordon term which is 
difficult to analyze. Therefore, some insight into the rel- 
evance of Hi s is helpful from a theoretical point of view. 
Following the philosophy of the previous subsection, we 
restrict the discussion to the spin polarization proper- 
ties of the ground state. The analysis of more compli- 
cated observables such as spin-spin correlation functions 
is beyond the scope of this work and will be discussed in 
another paper. 



Lowest eigen energies in different S z subspaces 
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FIG. 7. (Color online) The lowest eigen energies of different 
Sz subspaces with (solid blue) and without (dashed green) 
backscattering calculated for an edge of length L = 48. The 
broken SU(2) symmetry in the case of Ab s 7^ 1 lifts the de- 
generacy of the ground states in the different S z subspaces. 

Fig. [JJ compares the lowest eigen energies of the S z 
subspaces from calculations with and without backscat- 
tering. The most striking feature of the suppression of 
backscattering is the lifting of the ground state degener- 
acy in the SEM regime. This effect is easily understood 
by noting that the very reason for the ground state de- 
generacy in the Abs = 1 case was the SU(2) symmetry, 
which, however, is broken for Abs < 1- Interestingly, the 
lifting of the degeneracy is such that the lowest energy 
states of the subspaces with highest \S Z \, in the ground 
space for Ab s = 1, form the ground state for A bs < 1- 
This means that suppressing backscattering introduces 
an Ising anisotropy along the spin quantization axis cho- 



sen in the definition of the model. 

Apart from this degeneracy lifting, the evolution of 
the ground state properties with vpir/U is very similar 
for calculations with and without backscattering. The 
positions of the highest spin-polarization steps are prac- 
tically unchanged. Only at the steps close to the phase 
transition between the LL and the WEM regime, a de- 
viation of the Abs = results from the Ab s = 1 results 
can be observed. A handwaving explanation of this be- 
havior can be given in terms of the bosonization analysis 
of the WEM regime in Ref. [l3j. As soon as the Fermi 
levels for the up-spin electrons and the down-spin elec- 
trons are split, the backscattering process for electrons 
right at the Fermi surface is forbidden because it is not 
momentum-conserving. Thus, in order to conserve mo- 
mentum, the electrons are forced to scatter to higher en- 
ergies if there is a non-zero spin-polarization. This mech- 
anism suppresses the backscattering. In the bosonization 
language the backscattering Hamiltonian acquires a spa- 
tially oscillating phase which makes the corresponding 
operator irrelevant in the renormalization group. Thus 
in the WEM regime, not too close to the critical point, 
H^ s is suppressed and docs not give an important con- 
tribution. 

Close to the critical point, however, Fig. [7] indicates 
that Hi s becomes more important. This observation is 
consistent with the qualitative bosonization argument: 
Close to the critical point the phase oscillations in the 
bosonic backscattering Hamiltonian get slower until they 
completely disappear at the critical point. 



IV. DISCUSSION 

On the basis of a generalized class of effective models 
for one-dimensional interacting electrons we have stud- 
ied the magnetic properties of a graphene zigzag edge. 
Using exact diagonalization we confirmed the existence 
of three phases within these models, namely the satu- 
rated edge magnetism phase which is present at normal 
graphene edges, the Luttinger liquid phase which appears 
for edge states with strongly enhanced bandwidth, and 
an intermediate regime of weak edge magnetism. The 
latter phase is a realization of a ferromagnetic Luttinger 
liquid, a one-dimensional itinerant ferromagnet. We pre- 
sented evidence that the transition between the Luttinger 
liquid and the weak edge magnetism phase becomes a sec- 
ond order quantum phase transition in the limit of long 
edges. 

Beyond the identification of the magnetic properties 
of edge states, we examined the question why electrons 
in one-dimensional edge states have such a rich phase 
diagram with two types of ferromagnetic ground states, 
while usual one-dimensional electrons do not show any 
ferromagnetism. In view of the Lieb-Mattis theorem,— 
which actually forbids a spin-polarized ground state for 
interacting electrons in one dimension, this question be- 
comes even more pressing. 
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A closer inspection of the edge state model, which was 
derived directly from the graphene crystal structure^ 
revealed two unusual features of edge states which can- 
not be found in other one-dimensional electronic sys- 
tems, such as quantum wires, for instance. These are 
(a) the total absence of umklapp processes in the effec- 
tive electron-electron interaction, independently of the 
filling factor, and (b) a strong momentum dependence 
of the effective interaction vertex. Each of these fea- 
tures precludes the applicability of the Lieb-Mattis the- 
orem. The momentum-dependence gives rise to a com- 
plicated non-local interaction, which cannot be written 
as V(x, y)n(x)h(y), with h(x) an electron density oper- 
ator, as it is required for the Lieb-Mattis theorem. And 
even in the limit of momentum-independent interactions 
(ri = in the general model) the suppressed umklapp 
scattering makes the reformulation as a density-density 
interaction in real space impossible. In order to further 
track down the particular consequences of these special 
features for the magnetic properties of graphene edges, 
we replenished the direct model with a tunable umklapp 
scattering term and replaced the interaction vertex by 
a generalized vertex function in which the momentum- 
dependence can be switched on and off. 

The study of this generalized model, which can be 
tuned continuously between its edge state limit and a 
regime in which it describes normal one-dimensional met- 
als, revealed the significance of the two edge state fea- 
tures: the absence of umklapp processes is responsible 
for the existence of a spin-polarized ground state, and the 
strong momentum dependence of the interaction vertex 
stabilizes a regime of weak edge magnetism and gives rise 
to a second order phase transition between the param- 
agnetic Luttinger liquid and the ferromagnetic Luttinger 
liquid. 

It is interesting to note that the stabilization of the 
weak edge magnetism phase seems to be very robust 
against changes in the details of the interaction vertex 
function. Apparently it is only important that the ver- 
tex is suppressed as one of the four momenta of the par- 
ticipating fcrmions gets close to one of the Dirac points. 
The exact functional form of this suppression, however, 
seems to be irrelevant, since the qualitative behavior of 
the spin-polarization did not depend on whether we used 
the interaction vertex of the direct model or the interac- 
tion vertex of the generalized model with maximal mo- 
mentum dependence. These two vertex functions have 
in common that they vanish as one of the momenta ap- 
proaches a Dirac point. However, their functional forms 
are very different. 

Finally we note that one-dimensional itinerant mag- 
netism has also been studied in Hubbard chains with an 
additional second neighbor hopping , 30 ! 31 showing that 
it is indeed possible to define one-dimensional models 
which, at first sight, seem to comply with the prerequi- 
sites of the Lieb-Mattis theorem, but nevertheless have a 
high spin ground state. T he p hysical picture behind the 
model discussed in Refs. [3Q and Hi], however, is much 



different from the present work. Interestingly, the sign 
of the hopping amplitude to the nearest neighbor must 
be different from the sign of the next-nearest neighbor 
hopping for the system to have a ferromagnetic ground 
state. 

Also, we emphasize that the model discussed here 
is the low-energy theory of a realistic system which 
may be studied experimentally. It has been derived 
in direct line from a two-dimensional lattice model of 
graphenc/graphane interfaces ! 13 ' 20 
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Appendix A: Variational analysis of the generalized 
model 



We calculate the magnetic ground state properties 
of the generalized model within a fermionic mean-field 
approximation. It is assumed that only the averages 



k ra Ckro J are non-zero, so that the umklapp Hamilto- 

nian ff" s and the backscattering Hamiltonian H\ s drop 
out of the mean- field treatment. The resulting non- 
interacting Hamiltonian is diagonal in the momentum 
k. in the direction of motion r and in the z-spin pro- 
jection, so that the mean-field theory is equivalent to a 
variational ansatz based on the trial wave function 



m = n 



n 

.k<.kfT r 



-kRa 



n 

k~> — kp 



^kLo 



|0> (Al) 



with an asymmetric occupation of spin-up and spin-down 
states. The variational parameter M £ [0, 1] is related to 
the spin-dependent Fermi levels by 



^Fa =a-M, 
6 



(A2) 



and to the spin-polarization S, used in Sect. IIII1 by 
M = S'/S'max- For finite size systems, as discussed in 
the main part of this paper, the Fermi level cannot be 
varied continuously so that also M is a discrete variable 
in this case. However, within mean-field theory it is easy 
to perform the calculations in the thermodynamic limit, 
so that we will consider M to be a continuous variable 
and interpret it as the magnetization order parameter. 

The variational energy E(M) is easily calculated from 
the Hamiltonian H in Eq. (fTTj) 



E(M) = (M\ H \M) = -^(kvf - U)M 2 



5184 



(A3) 
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For U < tvv f , the minimum of E(M) is at M = 0, while 
for U > irvp the mean-field ground state has a finite 
magnetization 



M 



VpTX 



, 1 



(A4) 



Note that by definition the magnetization cannot become 
larger than 1. From Eq. (|A4[) it becomes obvious that 
a non-zero momentum-dependence Ti is required to sta- 
bilize the regime of weak edge magnetism. For Ti = 0, 
the magnetization would jump from to 1 at the critical 
point U = vpTT. 

The existence of the weak edge magnetism can be 
traced back to the M 4 term in Eq. (|A3|) which is gener- 
ated by the momentum-dependence Ti of the interaction 
vertex. In dimensions higher than one, such M 4 terms 
emerge also from momentum-independent interactions or 
directly from the kinetic energy, so that at least on the 
mean-field level Ti > is required for the stabilization 
of weak ferromagnetism only in one dimension. 



Appendix B: Exact diagonalization of the direct 
model 



The direct model Hamiltonian H dal defined by Eqs. 
(0 - 01 and the edge state limit of the generalized model 
Hamiltonian H [Eq. (fTTj) with Ab s = 1, ^us = and 
Ti = 6/tt] are not exactly equal, as the general model 
linearizes the single particle dispersion and replaces the 
factors Af p by the approximation SJ[. Nevertheless, the 
most important properties of graphene edge states, i.e. 
the momentum-dependence of the interaction vertex and 
the absence of umklapp scattering, are properly described 
by both, the direct model and the general model in the 
edge state limit. 

In this appendix, we check that the direct model has 
qualitatively the same magnetic properties as the general 
model in the edge state limit. In Fig. [8j we present 
the spin-polarization 5 as a function of A/U, obtained 
from the exact diagonalization of the direct model. The 
bandwidth parameter A of the direct model corresponds 
to the Fermi velocity vf of the general model. 



Clearly, for larger A/U, which corresponds to the pa- 
rameter vftt/U in the generalized model, we obtain a 
Luttinger liquid phase with a ground state of total spin 
S = 0. An intermediate regime with weak edge mag- 
netism exists, where the total spin of the ground state 
S < S'max is not maximal. As in the exact diagonalization 
analysis of the general model in the main text, only some 
of the lowest eigen energies in different S z subspaces are 
degenerate and form the ground state. For small A/17, 
the saturated edge magnetism phase is reached and the 
spin of the ground state is maximal, i.e., the lowest eigen 
energies in all S z subspaces are degenerate. 

Note that the energy of the fully spin polarized eigen- 
state of the direct model Hamiltonian has a finite slope 

Lowest eigen energies in different S z subspaces 
for the direct model 
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FIG. 8. (Color online) Lowest eigen energies in the S z = 
0, ±2, ±4, ±6, ±8 subspaces (from bottom to top) for the di- 
rect model calculated for an edge of length L = 48. The inset 
shows the dependence of the spin-polarization S as a function 
of A/U determined from the degeneracy of the ground state. 



(sec Fig. This is because the direct model lacks 

a symmetry of Hq of the generalized model leading to 
i?Q max (t;F) = const (cf. Fig. As only the degener- 
acy of the lowest eigen energies are important, but not 
their absolute values, this difference does not have any 
physical consequences. 
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